KANAZAWA-05-05 
RBRC-489 



The role of the double pole in lattice QCD 
with mixed actions 



Maarten Golterman," Taku Izubuchi 6 and Yigal Shamir c 

a Department of Physics and Astronomy, San Francisco State University 
San Francisco, CA 94132, USA 
maarten@stars . sfsu.edu 

b RIKEN-BNL Research Center, Brookhaven National Laboratory 
Upton, NY 11973, USA, and 
Institute for Theoretical Physics, Kanazawa University 
Kakuma, Kanazawa, 920-1192, Japan 
izubuchi@quark . phy . bnl . gov 

c School of Physics and Astronomy 
Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel-Aviv University, Ramat Aviv, 69978 ISRAEL 
shamir@post . tau .ac.il 

ABSTRACT 

We investigate effects resulting from the use of different discretizations 
for the valence and the sea quarks in lattice QCD, considering Wilson 
and/or Ginsparg- Wilson fermions. We assume that such effects appear 
through scaling violations that can be studied using effective lagrangian 
techniques. We show that a double pole is present in flavor-neutral Gold- 
stone meson propagators, even if the charged Goldstone mesons made out 
of valence quarks and those made out of sea quarks have equal masses. 
We then consider some observables known to be anomalously sensitive to 
the presence of a double pole. For these observables, we find that the 
double-pole enhanced scaling violations may turn out to be rather small 
in practice. 



1 Introduction 



Advanced numerical computations in Lattice QCD have recently seen an increase 
in the use of mixed actions. With "mixed action" we refer to any lattice theory in 
which differently discretized fermions are used in the valence- and sea-quark sectors. 
This development originates from the desire to use a relatively affordable method to 
simulate the sea quarks, using for instance various kinds of Wilson-like or staggered 
fermions, combined with a need for very good chiral symmetry in the valence sector, 
which can be obtained by using staggered or Ginsparg-Wilson quarks. 

The worry about mixed actions is that they do not satisfy some of the basic prop- 
erties of a healthy quantum-field theory, such as the possibility to analytically continue 
to Minkowski space and unitarity. We recall that the same is true for quenched and 
partially quenched (PQ) theories. However, the PQ theory reduces to full QCD if 
valence- and sea-quark masses are taken equal. In the mixed-action case, no such 
limit exists at non-zero lattice spacing. 

In this paper we will make the following working assumptions. (1) Any new 
unphysical effects resulting from the use of a mixed action disappear in the continuum 
limit; (2) At finite lattice spacing, these effects are controlled by positive powers of the 
lattice spacing a, i.e. they only show up through systematically manageable scaling 
violations; (3) One can use effective-field theory techniques to investigate their size. 

Effective field theory has proven to be a very useful tool for investigating the 
unphysical effects resulting from quenching and partial quenching, and we expect the 
same for any new effects resulting from the use of mixed actions. The idea is then 
simply to look for observables which are expected to be maximally affected by the 
scaling violations of a mixed action. Our experience with PQ theories gives us a clue 
for where to look. When the valence-quark masses are chosen unequal to the sea- 
quark masses, the most infrared-sensitive unphysical effects originate in the double 
pole that occurs in flavor-neutral Goldstone-meson two-point functions jl]. 

More generally, the double pole is an unphysical effect which arises from choosing 
the valence-quark action unequal to the sea-quark action. Other unphysical effects 
are known to arise from lattice-fermion actions containing beyond nearest-neighbor 
couplings, and from gauge actions containing Wilson loops larger than a single pla- 
quette. These actions often lack reflection positivity, and, hence, unitarity at finite 
lattice spacing. Under such circumstances it is imperative to take the continuum limit 
before analytically continuing to Minkowski space. It is generally believed that the 
notion of universality applies to the continuum limit obtained this way: one expects 
the same continuum limit regardless of the details of the regularization, i.e. of the 
lattice discretization. We expect that universality extends to mixed actions, and our 
working assumptions translate this expectation into more concrete terms. A practi- 
cal test then lies in comparing the results we obtain using chiral perturbation theory 
to numerical lattice calculations. Agreement would indicate that the assumptions 
underlying our effective-field theory calculations are valid. 

We will consider three classes of mixed actions, namely those with Wilson-like 
(Wilson, clover, twisted-mass Wilson) sea and valence quarks, those with Wilson-like 
sea quarks and Ginsparg-Wilson (GW) valence quarks, and those with GW valence 



2 



and sea quarks. The latter refers to any lattice fermion satisfying the GW relation [2] 
either exactly or to a good approximation; overlap, domain-wall, and perfect-action 
fermions included. 

In all cases, we will assume that different discretizations have been chosen in the 
valence and sea sectors, i.e. that the lattice Dirac operators for the valence- and the 
sea-quarks differ not only in their bare mass parameters (this includes for instance 
smearing the links in the valence sector but not in the sea sector). For simplicity we 
will restrict the discussion to degenerate masses in the sea and valence sectors. Our 
results can be generalized straightforwardly to non-degenerate and twisted quark- 
mass matrices. 

The chiral lagrangian for mixed-action lattice QCD with Wilson-like and GW 
fermions has been considered before in Refs. [SI E], which we take as our starting 
point. 1 In Sec. 2 we discuss the relevant elements of the chiral effective theory. We 
calculate the residue of the double pole in the two-point function of flavor-neutral 
Goldstone mesons, including the leading-order scaling violations. We then derive our 
main result: In a mixed-action theory there is a double pole in this two-point function, 
even if the charged Goldstone mesons made out of valence quarks and those made out 
of sea quarks have equal masses. If the sea and/or the valence action is Wilson-like, 
the residue of the double pole is of order a 2 . 

We then go on in Sec. 3 to consider the size of this effect for some interesting 
physical quantities which are particularly sensitive to the anomalous infrared behavior 
originating from a double pole. In particular, we consider the 1 = two-pion energy 
shift in finite volume, where the double pole causes an anomalous enhancement as a 
function of the volume [Oj, as well as the two-point function of the 1=1 scalar meson 
(the do), which has anomalous behavior for large euclidean time [7]. 2 Our conclusions 
are in Sec. 4. 

2 The residue of the double pole in mixed theories 

Our analysis starts from the results of Refs. [3 [I]. We first consider the implemen- 
tation of both valence and sea quarks through Wilson-like fermions. "Wilson-like" 
stands for any lattice fermion with exact flavor symmetry, 3 but with no chiral symme- 
try at non-zero lattice spacing a. In Ref. [1] first the Symanzik effective quark action 
was constructed 4 for Wilson fermions to order a 2 . If both valence and sea quarks are 
made of the same type of Wilson-like fermion, the (vector) symmetry group is that 
of a partially-quenched (PQ) theory pQ, SU(K + N\K), where K is the number of 
valence quarks, and N is the number of sea quarks. In the continuum limit chiral 

1 For very recent work considering GW valence quarks with a staggered sea, see Ref. [5]. Here we 
will not consider staggered fermions. 

2 Similar effects also occur in K° — > ir + ir~ matrix elements 0. 

3 As already mentioned, for simplicity we limit the discussion to degenerate masses in the sea and 
valence sectors. The generalization to non-degenerate and twisted |H] masses is straightforward. 
4 For a review of the method, as well as further references, see Ref. |10|. 
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symmetry is restored (for massless quarks), and the symmetry group enlarges to 5 

SU(K + N\K) L x SU(K + N\K) R . (partially quenched) (1) 

Ref. j3] also considered Wilson-like sea quarks with valence quarks of the GW 
type E] ■ The a/0 mixed-action theory is then invariant under 

[SU(K\K) L x SU{K\K) R ] x SU(N) . (GW - Wilson) (2) 

Again the complete Symanzik effective action was constructed to order a 2 . For both 
cases, the complete chiral lagrangian was derived to the same order from the Symanzik 
effective action with the help of a spurion analysis. 

Here we are only interested in the role of the double pole in neutral Goldstone 
meson propagators, and we will restrict ourselves to operators in the chiral lagrangian 
relevant for our purpose. However, our aim is to consider also mixed actions in which 
both valence and sea quarks are Wilson-like, but different. In that case, the symmetry 
group is reduced to 

SU(K\K) x SU(N) , (Wilson - Wilson) (3) 

and more operators appear in the effective theory. 

In the continuum limit, the chiral effective lagrangian is just that of PQ QCD 
PP. In terms of the non-linear field £ = exp(2z$//), with <3> the matrix of Goldstone- 
meson fields, the lowest-order (euclidean) chiral lagrangian is 

£ C ont = \f str(0 M £0 M £t) - l - f 2 B str(£A^ + M^) + l - /i 2 $ 2 , (4) 

where <3> = str($) = (f/2i) logsdet(£), and /, B and /i 2 are low-energy constants 
(LECs). 6 M. is the quark mass matrix, 

M = M ] = di&g(m v , m v ,..., m s , m s , . . . , m v , m v , . . .) , (5) 

while str and sdet are the supertrace and the superdeterminant needed to construct 
invariants for graded groups. The first K masses are those of the valence quarks, the 
next N those of the sea quarks, and the final K those of the ghost quarks needed to 
cancel the loops of the valence quarks [T4] . 

In Eq. (0J) we kept the singlet mass term because we are interested in scaling 
violations to the singlet mass parameter /i 2 . Normally one would not do this because 
fj, is considered as being of same order as the chiral-lagrangian cutoff scale A w 47r/. 
The pseudo-scalar singlet's mass-squared is equal to Nfi 2 + 2Bom s , and so this degree 
of freedom can be integrated out. Accordingly, one could constrain £ by setting 
sdet(£) = 1, i.e. $o — 0. We will indeed recover this result in what follows, and from 
then on limit ourselves to the theory with the pseudo-scalar singlet integrated out. 
This is also the reason that we did not keep any other terms containing the singlet 
field $ in Eq. (gj). 

5 For explanations about the legitimacy of using this symmetry group, see Refs. |11M12| . 
6 Note that /i 2 (this paper) = /i 2 (Ref. Q])/3. 
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Because of the reduced symmetry group (jUJ) for a Wilson-like mixed-action, the 
following operators appear at order a 2 in the chiral lagrangian: 



5£w = -^( 7 ™(str(P„(S-St))) 2 + 7ss (str(P s (S-St))) 2 

+2 lvs stx(P v (E - str(P s (£ - E+))) . (6) 

Here P v and P s are projectors on the valence, respectively, sea sectors, 

P v = diag(l„, 0, lg) , P s = diag(0, l s , 0) , (7) 

where l v and l g are the KxK identity matrices in respectively the valence and ghost 
sectors, and l s the N x N identity matrix in the sea sector. Eq. (jUJ) follows from an 
analysis very similar to that of Ref. It was shown there that the corresponding 
operators in the Symanzik effective action are four-fermion operators, which appear 
at order a 2 . As usual, "translation" to the chiral effective lagrangian is done using 
a spurion analysis. Here the spurions are needed to restore both chiral symmetry 
and full valence-sea symmetry. When expectation values are assigned to them, both 
these symmetries are explicitly broken, cf. Eq. (JHJ). 

In Ref. |3] only one similar operator appears in the chiral lagrangian, because the 
same Wilson-fermion action was assumed for both the valence and sea quarks. The 
symmetry group is then SU(K + N\K) instead of (J3J), which implies 7 

1w = Jvs = Iss ■ (partially quenched, Wilson) (8) 

We will also consider mixed actions in which the valence quarks are of the GW 
type, while the sea quarks are Wilson-like [3]. The exact chiral symmetry in the 
valence sector (cf. Eq. (J2J) implies j vv = r ) vs = in Eq. (JBJ). What survives is only 
the term corresponding entirely to the sea sector, 8 where chiral symmetry is still 
broken at a ^ 0. 

Finally, when both valence and sea quarks are of the GW type, no such operators 
can appear at all, and all three 7's vanish, consistently with the enlarged symmetry 
group (for massless quarks) 

[SU(K\K) L x SU(K\K) R \ x [SU(N) L x SU(N) R \ . (GW - GW) (9) 

Operators similar to those in Eq. (jUJ) can appear only at order a 2 m 2 , with m ~ m v ~ 
77V At the level of the chiral lagrangian, they can be made by replacing the factors 
str(P WiS (S - St)) in Eq. © by str(P 1 , jS (SA^ t - Mtf)), where M is now promoted to 
a spurion field transforming like S. 

For domain- wall fermions with a finite fifth dimension L5, chiral symmetry is only 
approximate (A mixed action could result e.g. from choosing different values for 
L5 in the valence and sea sectors.) The ensuing chiral-symmetry violating effects can 

7 This corresponds to the operator with LEC W!j in Eq. (21) of Ref. 4 . The relevant spurions 
are B and C of Eq. (14) of Ref. 4 . Replacing the identity matrix I in that equation by P v or P s in 
all possible combinations leads to our Eq. 

8 The corresponding operator is the one multiplying Wy in Eq. (40) of Ref. 0]. 
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be taken into account by adding separate spurions for the residual quark masses 
uni nzi- As long clS TYl r Q S IS smaller than the explicit quark masses, 9 again all such 
operators are of order a 2 m 2 . 

In order to obtain tree-level meson propagators, we expand £ ccm t + ££w t° 
quadratic order in the meson field $. It is useful to define separate singlet fields 
in the valence and sea sectors: the "super-?]'" $q of the valence-symmetry group 
SU(K\K) [IH1, and the "dynamical" singlet field T]' s . Explicitly, 

$0= £ £ *u, v' s = £*«, ( 10 ) 

i valence i ghost i sea 

where $o = $0 + V's- The singlet mass terms following from Eqs. ()4I6|) can now be 
written as ^ ^ 

^singlet = ^ /4 (*o) 2 + ti. %V's+2 »ls ^'sf > ( U ) 

in which 

l4v = I? + IwO 2 , (12) 
l4s = V 2 + Ivsa 2 , 
ffis = ^ + l ss a 2 . 

These scaling violations contain no derivatives, nor (for Wilson fermions) powers 
of the quark mass matrix. Hence they are part of the leading-order difference in 
the chiral expansion between a mixed action at non-zero lattice spacing and its PQ 
continuum limit. 

Before proceeding, we have a few comments. First, new operators as in Eq. © 
lead not only to order a 2 contributions to singlet mass terms, but also to new 0(a 2 ) 
vertices in the chiral lagrangian. In other words, chiral symmetry implies that Eq. (fTTj) 
cannot appear as an isolated invariant term. 10 Second, while other operators appear 
at order a 2 in the chiral lagrangian |E|, all the quadratic terms are of the form 

a 2 (p v str( J R u $ 2 ) + f3 s str(P s $ 2 )) , (13) 

with /3 V)S linear combinations of order a 2 LECs. 11 They do not contribute any singlet 
mass terms, i.e. terms as in Eq. (fTTj) . 

With our quark mass matrix (JHJ), all pions made out of valence quarks are de- 
generate with a common mass M vv , and, similarly, pions made out of sea quarks all 
have a common mass M ss with, in general, M ss 7^ M vv . Here M vv denotes the mass 
of a charged valence meson, i.e. a meson made out of a quark and an anti-quark of 

9 If L5 is chosen so small {e.g. L5 = 2) that the residual masses, and, hence, the actual quark 
masses are substantially bigger than the explicit quark-mass parameters, domain-wall fermions 
should be treated as Wilson-likc. 

10 We thank Claude Bernard for very useful discussions on this point. 

n In the GW-valence, Wilson-sea case, j3 s is a linear combination of Wg and W% in Eq. (40) of 
Ref. gj. 
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different ffavors. At tree level, these meson masses are given by [4 12 

M 2 V = 2B 0v m v + 2W 0v a + 2P v a 2 + 0(am,m 2 ,a 3 ) , 
Ml = 2B 0s m s + 2W 0s a + 2f3 s a 2 + (am, m 2 , a 3 ) , 



where Bq v — Bq s = 0(a 2 ). The LECs Wq v or Wq s vanish if one works with an order-a 
improved action in the valence or sea sector. The terms of order a and a 2 both vanish 
if one uses fermions of the GW type. A key observation is that, at tree level, the 
operators in 5£\y only contribute to singlet mass terms, as in Eq. (JTTJ). The LECs 
1w,vs,ss do not appear in the tree-level expressions for the meson masses M vv>ss . 

We obtain the tree-level two-point function Gij = $^-) for the neutral valence 
meson fields as 

r n &ij vL (P 2 + M i) + N{fi 2 vvf i 2 s - fits) n ,x 
hAP) P 2 + M 2 V (p' + MlW + Ml + Nvl) ' 1 1 

where the indices i and j correspond to the flavor of the valence quark. Since the 
effective theory is based on a combined expansion in M 2 V ss ~ p 2 and a 2 , treating fi 2 
as a quantity of order one in this expansion, we can simplify this expression to leading 
order in these parameters by using Eq. (J12|) . We obtain 13 

GlM = {*« ~ n) WTmI {p 2 + miy ' (16) 

Indeed, /i 2 disappears from the propagator, consistent with the fact that in a PQ 
theory $ = str($) can be constrained to vanish [20] • Setting $ = implies that 
r]' s = an d t nus £ S i n giet goes over into 

Aingiot ^ fl2 (7™ + 7ss - 27„ s )($o) 2 • ( 17 ) 

This simplified singlet mass term directly leads to Eq. (|16jl. Setting str($) = also 
removes terms of the form $ str(P„ )S (E — H*)), ((9 M $o) 2 , etc., so that we were justified 
in not considering such terms in Eq. (pj) . 

From Eq. ()16|) . we learn that in a mixed-action theory, double poles may arise from 
two different sources. One is well known from PQ QCD. It has a residue proportional 
to the difference between valence- and sea-meson's mass-squared. In addition, at non- 
zero lattice spacing, a new term proportional to a 2 appears in the residue. This new 
term arises because scaling violations in the valence, sea and mixed double hairpins 
are different for such theories, and thus 7^^ + 7 SS — 2 , y vs will in general not vanish. 
Of course, if we use the same lattice fermions for both valence and sea sectors, we 
recover the PQ case in which this new contribution to the residue does vanish, cf. 
Eq. ©. 

12 We assume that 1/a divergences have been absorbed into m„ iS . We also assume that we are 
close to, but outside the Aoki phase |19j . 

13 Note that we do not specify the relative size of M% v ss and a 2 . 
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To the order we work, it is possible to tune the residue of the double pole to 
zero by adjusting the valence quark mass. A vanishing residue would be obtained by 
choosing the valence quark mass m v such that the residue, 

R = (Ml - M 2 VV )/N + a 2 ( lvv + lss - 2 lvs ) , (18) 

vanishes. Alternatively, one could tune the valence quark mass such that M vv = M ss . 
At non-zero lattice spacing, these two choices are different. In other words, while also 
the meson masses are afflicted by scaling violations (cf. Eq. (|14jl). there are scaling 
violations in the double pole which do not get absorbed into the meson masses. If one 
chooses to tune the valence quark mass such as to set R — 0, obviously this requires 
a quantity sensitive to the singlet part of the rj' mass-squared, N if 2 . 

Our analysis does not apply to any mixed action involving staggered fermions. 
Mixed actions with only staggered or GW-type fermions have no order-a scaling 
violations. Order-a 2 scaling violations associated with symmetry breaking can only 
come from the breaking of taste symmetry in the staggered-fermion sector. As a 
result, all the order-a 2 scaling violations in the double pole can be related to the 
order-a 2 mass splittings between the staggered pseudo-scalar mesons For mixed 
actions with staggered and/or GW-type fermions, scaling violations that cannot be 
absorbed into the Goldstone-meson masses occur in the residue of the double pole 
only at order a 2 m 2 . 

3 Examples of quantities with enhanced scaling vi- 
olations 

The occurrence of a double pole in neutral-meson propagators leads to an enhanced 
sensitivity of many quantities to the meson masses [TBI ED HI E2| or to the volume 
[SI EJ EHj • This implies ipso facto an enhanced sensitivity to scaling violations of such 
quantities because of the presence of the LECs "im^s^s- In this section we discuss 
examples of such quantities, in particular the 1 = pion scattering length j£l and the 
two-point function for the 1 = 1 scalar meson a E2] • 

The oldest examples of enhanced sensitivity due to the double pole are the en- 
hanced ( "quenched" ) chiral logarithms in Goldstone meson masses, decay constants, 
and the chiral condensate fHl EJJ ■ However, for quark masses used in typical simula- 
tions, these effects are numerically very small; the same is expected to be true for the 
scaling violations present in Eq. ()16|) . A more interesting quantity is the shift in the 
energy of two pions, both at rest in an / = state, in a finite spatial volume L 3 with 
periodic boundary conditions. In quantities like this, the presence of a double pole 
gives rise to enhancement of finite- volume effects jSj. It is therefore a priori not clear 
what the combined effect is of enhancement by powers of L in addition to suppression 
by powers of a. 

The two-pion 1 = energy shift is given by jB] 

W = + 5*<*- t >* +1 2 i 1 - 7f) A "^ L ^ + °^ ' < 19 > 
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in which 

a = _ = M ' vv or\\ 

~ 87T 2f2 > " 167r 2 /2 > 

and R is defined in Eq. (|18|1. The first term is the leading-order, tree-level result, 
while the other terms come from one-loop diagrams. We adapted the result of Ref. [0] 
to the PQ case, resulting in the extra —1/N appearing in the Aq term. 14 

The large- L expansions of the coefficient functions Bq and Aq are given by 15 



B„(ML) = -^ + ^(^1) +0((£L) ], (21) 




5z(3) 5z(2) 
8tt 2 + 32^ 

° V ; 4vr 2 

in which z(3) and z(2) are geometrical constants, 

z(3) = 8.40192397 , z{2) = 16.53231596 . (22) 

We only showed the terms that generate enhanced finite-volume corrections to the 
leading-order result. These corrections are unphysical, and are absent if exactly the 
same lattice-fermion action is used in the valence and sea sectors (with valence-quark 
masses equal to sea-quark masses). When a mixed action is used, if one does not 
tune m v to make the residue in Eq. (fT%|) vanish, AEj =0 would be afflicted by scaling 
violations. (As already explained, this is true even if m v is tuned to set M vv = M ss .) 
The enhanced scaling violations are of order M vv a 4 L 3 Aq CD (from the leading term in 
B ), respectively M vv a 2 LAq CD (from the leading term in A ), relative to the physical 
value (estimating / ~ Aqcd and 

r Yw,vs,ss ~ Aqcd)- There are additional four-point 
vertices of order a 2 which contribute to these one- loop terms in AEj = q, that were 
not taken into account in Ref. jUj . However, they will not lead to any parametrically 
larger contributions than the ones considered here. 

It is instructive to consider these results in view of various power-counting schemes 
proposed in the literature. Ref. [3] considered an expansion with a single small pa- 
rameter (compare Eq. (J2UJ)) e ~ M 2 VSS /A 2 ~ p 2 /A 2 ~ aA QCD . Assuming M VV L > 1 
and M vv = M ss , the first term of Eq. (jTUj) is of order e/(M vv L) 3 , the second term 
is of order e 4 , while the third term is of order e 3 /(M vv L) 2 . If we would ignore the 
enhancement of the second and third terms by positive powers of M VV L relative to 
the first term, these terms would be of higher order compared to the 0(e 2 ) term in 
Eq. (fTT^) and the 0(e 2 ) terms considered in Ref. j3j. (This is consistent with the fact 
that, with this power counting, the a 2 term in R is of higher order than the M 2 S — M 2 V 
term.) In a different scheme [23], (oAqcd) 2 is taken to be of order e. In that case, 
the second and third terms in Eq. ()19j) are both of order e 2 . Again they are enhanced 
(by (M VV L) 3 , respectively M VV L) relative to the other 0(e 2 ) terms in Eq. (JTHJ). 

For parameters used in realistic simulations, the volume-enhanced unphysical 
effects may nevertheless turn out to be rather small. For instance, setting aM vv ~ 0.2, 

1 In the fully quenched case, which corresponds to taking M^ s — > oo in Eq. (JTSJ, the corresponding 
factor is equal to one. 

15 This corrects a mistake in B (ML) in Eq. (18) of Ref. 0. 
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°Aqcd ~ 0.1 and L/a = 32 leads to scaling violations of about 0.7% (from B (M VV L)) 
and 6% (from A (M VV L)), relative to the first, physical term in Eq. (jTUj) . Note that 
the numerical coefficient of (M VV L)~ 2 in A (M VV L) is rather large, about 50. 

We see that, if Wilson-like fermions are used in the valence and/or sea sectors, the 
unphysical effects in AEj =0 are expected to be rather small in realistic simulations, 
although probably not entirely negligible. Of course, better knowledge of j vv + 7 SS — 
2 r y vs is needed in order to move beyond our crude estimates. 

As explained earlier, when GW fermions with exact chiral symmetry are employed 
in both the valence and sea sectors, the scaling violations (in R) are further suppressed 
by a factor of order (m/Aqco) 2 - Thus, the unphysical effects will be negligibly small 
if the quark masses are light enough. The same is true for domain-wall fermions, 
provided the residual mass is not larger than typical explicit quark masses (as easily 
accomplished in current simulations). 

Another quantity which may be strongly afflicted by the double pole is the 1 = 1 
scalar (ao) two-point function |7j. This effect was considered in the PQ theory |2~2~j . 
which we will generalize to the mixed-action case. We note that the sensitivity of the 
a two-point function was used to tune a mixed action theory to the full unquenched 
QCD limit in Ref. j2Sj- However, the scaling violations present in the double pole, 
Eq. (fTKj). were not considered in that study. 

The 1 = 1 scalar two-point correlator at vanishing spatial momentum was calcu- 
lated in the PQ theory [22J. Adapting that result to our case, the dominant contri- 
bution for large time t is given by 16 



C(t) 



m 

2L 3 



e -2M V3 t pr e -2M vv t i 

{(M 2 v + M 2 s )/N + RM vv t 



M 2 2 M 2 M 2 

vs vv vv 



(23) 



Again, R is the residue of the double pole (no particular assumption on the value 
of M 2 S - M 2 V is made). M vs is the mass of a meson made out of a valence and a 
sea quark, with M 2 S = (M 2 V + M 2 s )/2 to leading order. We have assumed that the 
mass of the r]' s , given by M 2 , = Nfi 2 + M 2 S , as well as the ao mass, are much heavier 
than any Goldstone-meson mass. For large enough t the linear term in t dominates. 
(Depending on the actual values of the M's and 7's in Eq. (fTSj) . C(t) could even 
become negative.) Setting M vv = M ss reduces Eq. (|2l3|) to 

B 2 e~ 2M ^ (N 2 ( lvv + lss -2 lvs )a 2 t \ 
C[t) 2-V^CT \~2~N K v ) ■ (24) 

The coefficient of the factor linear in t (measured in lattice units) is of order (oAqcd) 4 / 
(aM vv ), which is again quite small for a typical simulation. (In the special case N = 2 
the first two terms in Eq. (j2*lj) cancel each other, and the last term should be compared 
with the a Q and 7/ terms which are present in the full correlator.) 
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Our LEC B corresponds to 2/i of Ref. [221 7 an d our [i 2 corresponds to TOq/3. 
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4 Conclusions 



We have investigated two quantities which are expected to be particularly sensitive to 
the unphysical effects resulting from a mixed lattice-QCD action. We have assumed 
that there is no problem of principle, i.e. we have assumed that all such effects are 
encoded in scaling violations. Still, it is important to estimate their numerical size for 
realistic choices of MC simulation parameters. Using effective field theory techniques, 
building on the work of Refs. [3J Ej , we identified the double pole in the flavor- neutral 
Goldstone meson propagators as the most infrared- sensitive probe of these effects. 

We considered two examples of observables where the double-pole scaling viola- 
tions are parametrically enhanced, either by powers of M VV L (i.e. of the spatial size) 
in the case of I = two-pion energy shifts, or linearly with euclidean time t in the 
case of the ao propagator. We found that these scaling violations are likely to be 
rather small in currently realistic simulations. While this is good news for anyone 
interested in using mixed actions, care should be taken to monitor their size in prac- 
tice. Indeed, for the two-pion energy shift the effect — while small — is not negligible 
already according to our crude estimate, if a Wilson-like action is used in the sea 
and/or valence sector. 

The magnitude of these scaling violations depends on the observable in question. 
It could well be that other interesting observables exist for which the effect is larger. 
As another example, let us mention the nucleon-nucleon potential, where the contri- 
bution from one-pion exchange is given by a Yukawa potential. In PQ QCD with 
unequal valence- and sea-quark masses, the long-distance part of the potential gets 
an additional contribution from the double pole through rj exchange [21)]. 17 Including 
the scaling violations of Eq. (j!6|) leads to a potential of the form 

K W = „ w . v „ B . v ( 9 >M_M _ 9 ,|_) , (25) 

where we took N = 2. Here g& (go) is the / = 1 (/ = 0) axial coupling. The 
indices (1) and (2) refer to the two nucleons, the cr's to spin and the r's to isospin. 
The second, unphysical term clearly dominates the first term at large distances. The 
nucleon-77 coupling g is poorly known, but it could be larger than g^, thus decreasing 
the distance at which the unphysical term sets in. 

If lattice fermions satisfying the GW relation (or approximations thereof such as 
domain-wall fermions with finite fifth extent) are used in both the valence and sea 
sectors, chiral symmetry leads to an extra suppression of the unphysical mixed-action 
effects, by a factor of order (iti/Aqcd) 2 - 

We conclude with a practical comment. For Goldstone-meson masses in the few 
100 MeV range, one has M vv /A QC b ~ M SS /A QCB = 0(1). In comparison, (oAqcd) 2 = 
0(0.01). Therefore, when Wilson-like fermions are involved, one might have to control 
the difference M% v — M^ s at the few percent level in order to be able to single out 
the scaling-violations part of the residue R. This is likely to require the use of an 
observable with enhanced sensitivity to the double pole. If both sea and valence 

17 We thank Paulo Bedaque for reminding us of this example of quenching effects. 
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fermions are of the GW type and the quark masses are small enough, the effect could 
be too small to be detected. 
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